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Abstract. We use Aura/MLS stratospheric water vapour 
(H 2 O) measurements as tracer for dynamics and infer inter- 
annual variations in the speed of the Brewer-Dobson circu- 
lation (BDC) from 2004 to 201 1. We correlate one-year time 
series of H 2 O in the lower stratosphere at two subsequent 
pressure levels (68hPa, ~ 18. 8 km and 56hPa, ~ 19.9 km at 
the Equator) and determine the time lag for best correla- 
tion. The same calculation is made on the horizontal on the 
100 hPa (~ 16. 6 km) level by correlating the H 2 O time se- 
ries at the Equator with the ones at 40° N and 40° S. From 
these lag coefficients we derive the vertical and horizontal 
speeds of the BDC in the tropics and extra-tropics, respec- 
tively. We observe a clear interannual variability of the verti- 
cal and horizontal branch. The variability reflects signatures 
of the Quasi Biennial Oscillation (QBO). Our measurements 
confirm the QBO meridional circulation anomalies and show 
that the speed variations in the two branches of the BDC are 
out of phase and fairly well anti-correlated. Maximum as- 
cent rates are found during the QBO easterly phase. We also 
find that transport of H 2 O towards the Northern Hemisphere 
(NH) is on the average two times faster than to the South- 
ern Hemisphere (SH) with a mean speed of 1.15 ms -1 at 
100 hPa. Furthermore, the speed towards the NH shows much 
more interannual variability with an amplitude of about 21 % 
whilst the speed towards the SH varies by only 10 %. An am- 
plitude of 21 % is also observed in the variability of the as- 
cent rate at the Equator which is on the average 0.2 mm s _1 . 


1 Introduction 

The very low amounts of water vapour (H 2 O) measured 
in the lower stratosphere by Brewer (1949) implied that 
air masses must be entering the stratosphere through the 
cold tropical tropopause because nowhere else temperatures 
are low enough to allow for sufficient dehydration. Dobson 
(1956) used ground-based measurements of total ozone and 
observed minimum values in the photo-chemical source re- 
gion of the tropics and maximum values in the extra-tropics. 
Hence, he assumed a transport from the tropical stratosphere 
towards the pole in both hemispheres marked by low-latitude 
ascent and high-latitude descent. The assumption of a pole- 
ward residual circulation now called Brewer-Dobson circu- 
lation (BDC) is still valid today. Figure 1 shows the dif- 
ferent branches of the BDC (adapted from Plumb, 2002). 
Air masses ascend inside the tropical pipe, which represents 
a mixing barrier for tropical and extra-tropical air (Plumb, 
1996). The rising air is distributed symmetrically towards 
both hemispheres in the lowermost stratosphere (black ar- 
rows) and brings H 2 O from low- to high-latitudes. The lower 
stratospheric meridional component is usually called shallow 
branch of the BDC and is faster during summer (Bonisch 
et al., 2009) when air masses are quasi horizontally mixed 
out of the tropics. Mixing into the tropics is slow com- 
pared to mixing within mid-latitudes (Volk et al., 1996) high- 
lighted by a smaller arrow pointing back towards the Equator. 
In the upper stratosphere air masses are predominantly di- 
rected towards the winter pole where they mix and descend 
into the so called “surf zone” and the polar vortex. This is 
called the deep branch of the BDC which is relatively slow 


Published by Copernicus Publications on behalf of the European Geosciences Union. 


4564 


T. Flury et al.: Variability of the Brewer-Dobson circulation 



Fig. 1. Schematic of the Brewer-Dobson circulation in the strato- 
sphere. In the lowermost stratosphere equatorial air is transported 
and mixed towards both poles (thick black arrows) before it re- 
enters the troposphere at higher latitudes and returns towards the 
Equator (Eq). The rising air masses in the tropics are transported to- 
wards the winter pole higher up, cross the tropical pipe and descend 
into the “surf zone” and the polar vortex. We focus on the three 
branches highlighted by the thicker blue and black arrows which 
distribute H2O in the stratosphere. The ellipse and the two poleward 
directed arrows represent the Hadley circulation in the troposphere. 


(Plumb, 2002). Stratospheric air returns to the troposphere at 
higher latitudes and is ultimately brought back into the trop- 
ics (equatorward arrows inside troposphere), where it joins 
the tropospheric Hadley-cell with upwelling in the tropics 
and descent in the subtropics to form the subtropical ridge. 

Haynes et al. (1991) explained that the BDC is driven 
by breaking of high-latitude vertically propagating planetary 
and gravity waves which slow down the stratospheric zonal 
mean flow and induce a meridional residual circulation. The 
planetary wave breaking occurs mainly in the “surf zone” de- 
picted in Fig. 1. However, more and more evidence from re- 
analysis data and model studies suggests that the upwelling 
in the lower tropical stratosphere is a combination of mid- 
to high-latitude wave breaking and the dissipation of tropi- 
cal waves. Garny et al. (2011) show that slow and large scale 
upwelling is driven by wave forcing both in tropics and extra- 
tropics. Boehm and Lee (2003) state that TTL upwelling 
is generated as a response to poleward propagating Rossby 
waves that are driven by tropical convection. Randel et al. 
(2008) found that both equatorial planetary waves and waves 
originating from the extra-tropics create horizontal eddy mo- 
mentum flux convergence which forces upwelling across the 
tropical tropopause. Using a simple model Chen and Sun 
(2011) showed that the planetary wave forcing in the win- 
ter tropical and subtropical stratosphere contributes most to 
the seasonal cycle of tropical stratospheric upwelling rather 
than the high-latitude wave forcing. 


On the other hand Randel et al. (2002a, b), Ueyama 
and Wallace (2010) and Dhomse et al. (2008) showed di- 
rect evidence of the correlation of high-latitude planetary 
wave breaking and the strength of the BDC by using mid- 
latitude eddy heat flux as a proxy for wave breaking. A 
stronger BDC leads to lower tropical tropopause tempera- 
tures (Randel et al., 2006) and to stronger freeze drying of 
the stratosphere (Dhomse et al., 2008) i.e., lower H 2 O con- 
centrations. Knowledge of the BDC is, thus, crucial to deter- 
mine the budget of H 2 O and other atmospheric trace gases 
in the stratosphere since it has an important influence on the 
stratospheric entry, the residence time and ozone depletion. 
The BDC has been varying since observations started in the 
1960s (Roscoe, 2006). General circulation models (GCM) 
predict that in a warmer climate the mass exchange between 
the tropics and extra- tropics will increase (Butchart et al., 
2006; McLandress and Shepherd, 2009; Okamoto et al., 
2011). The predicted enhancement of the BDC is due to an 
increase in planetary wave drag in the extra-tropical strato- 
sphere as well as to an upward shift of the orographic gravity 
wave drag at low and mid-latitudes. A change in the BDC 
will in turn also affect stratospheric H 2 O and ozone abun- 
dance. 

Observational studies about the BDC variability are rare. 
Engel et al. (2009) showed by using balloon-borne measure- 
ments of stratospheric trace gases that the stratospheric mean 
age of air was slightly increasing rather than decreasing al- 
though no statistically significant trend could be extracted. 
This result disagrees with model predictions for a strength- 
ening of the BDC and an associated decrease in the mean age 
of air. Several previous studies used the H 2 O tape recorder 
signal to derive the ascent rate in the tropical lower strato- 
sphere as a measure for the upwelling branch of the BDC 
(Mote et al., 1996; Schoeberl et al., 2008; Niwano et al., 
2003; Fujiwara et al., 2010). The tropical H 2 O tape recorder 
is displayed in Fig. 2 and shows the imprinted seasonal cycle 
of H 2 O around the 100 hPa level which is then transported 
upward by the slow BDC ascent. Fujiwara et al. (2010) and 
Niwano et al. (2003) showed the seasonality and a clear 
modulation of the ascent rate by the Quasi Biennial Oscil- 
lation (QBO) using radiosonde measurements and HALOE 
satellite measurements, respectively. The QBO is a dynamic 
phenomenon in the equatorial stratosphere and manifests it- 
self by a quasi periodic oscillation (=28-29 months) of the 
zonal wind between westerlies and easterlies. The alternat- 
ing winds propagate downward from the stratopause (1 hPa, 
~55 km) at a rate of about 1 km per month towards the trop- 
ical tropopause (100 hPa, ~ 16.6 km). The ascent rate of the 
tape recorder is modulated by the QBO and speeds up dur- 
ing the colder easterly phase and slows down during the 
warmer westerly phase (Niwano et al., 2003; Plumb and Bell, 
1982). The ascent rate is thus anti-correlated with the temper- 
ature. As a consequence the strength of the BDC determines 
stratospheric H 2 O entry values by modulating TTL temper- 
atures. The QBO affects also the meridional wind. In fact, it 
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Fig. 2. The so called atmospheric tape recorder shows zonal mean 
H 2 O anomalies in the tropics as a function of time and altitude. 
Red colours show above average volume mixing ratios, blue shows 
lower values mainly during NH winter near lOOhPa. The seasonal 
cycle with an amplitude of about 30 % is imprinted at lOOhPa and 
slowly transported upward by the ascending branch of the BDC as 
indicated by the arrow. The slope of the arrow indicates the average 
speed of the tropical up welling. The variability between 316 and 
150 hPa is mostly due to ENSO. In El Nino (La Nina) years more 
(less) H 2 O is observed. 



induces a secondary meridional circulation due to thermal 
wind balance, which affects the meridional branch of the 
BDC (Plumb and Bell, 1982; Baldwin et al., 2001; Ribera 
et al., 2004; Punge et al., 2009). More about the QBO can be 
found in Baldwin et al. (2001) and references therein. 

The above mentioned work focused on seasonal and in- 
terannual variations in the tropical upwelling branch of the 
Brewer-Dobson circulation. By using Aura/MLS measure- 
ments we add observations of the speed of the shallow merid- 
ional branch at lOOhPa and concentrate on its link to the 
tropical upwelling branch of the BDC as highlighted by the 
thick arrows in Fig. 1. We use the rate of horizontal mixing 
and transport of H2O to mid-latitudes as a proxy for the shal- 
low meridional branch. For simplification we will call it the 
meridional speed of the BDC although it has to be pointed 
out that the speed derived from the meridional propagation of 
H2O anomalies does not provide an absolute measure for the 
net transport of mass by the BDC. This is in contrast to the 
tape recorder signal which consists primarily of mass advec- 
tion. The horizontally propagating anomalies we observe are 
due to a combination of residual mass transport and two-way 
mixing. Hence, the derived speed depends on the used tracer, 
but is a qualitative measure of the BDC. However, poleward 
residual mass transport is dominating. Trepte et al. (1993) 
used SAGE II satellite observations of stratospheric aerosols 
from the Mount Pinatubo (15° N, 120° E) eruption in 1991 to 


characterise the stratospheric circulation. Their results show 
that aerosols were predominantly moved poleward and show 
no signs of mixing back to their tropical origin in the lower 
stratosphere. 

Flury et al. (2012) showed that total water (ice and vapor) 
is roughly constant in the tropical tropopause layer (TTL) 
at lOOhPa on seasonal time scales. However, the tempera- 
ture determines the balance between ice in cirrus clouds and 
water vapour. Thus, the strong seasonal cycle of temperature 
manifests itself also in water vapour and ice and can be found 
in water vapour throughout the lower stratosphere from trop- 
ics to mid-latitudes due to transport by the BDC. The slow 
residual transport of air out of the lower tropical stratosphere 
leads to a time lag in the H2O time series at higher lati- 
tudes and altitudes because lower stratospheric H2O concen- 
trations are set in the tropics. This time lag enables us to re- 
trieve the average speed of the BDC thanks to the tracer char- 
acteristics of H2O. The method allows us to investigate the 
difference between the Northern (NH) and Southern Hemi- 
sphere (SH) by considering the different amplitude of the 
interannual variability of the respective meridional branch. 
Furthermore, we link the derived speeds to temperature and 
zonal wind variability. The article is organised as follows. 
Section 2 describes the utilised satellite and reanalysis data as 
well as the correlation method. Section 3 outlines the results, 
especially the interannual speed variations and their link to 
temperature and zonal wind. Section 4 discusses the QBO’s 
effect on the tropical pipe and how we interpret the results. 
Finally we conclude and summarise our findings in Sect. 5. 


2 Data and method 

In this study, we analyse measurements from the MLS (Mi- 
crowave Limb Sounder) instrument which is flown on board 
the Aura spacecraft and part of the NASA A-train satellite 
constellation. MLS measurements began in August 2004 and 
continue since. We use daily mean H2O, ozone (O3) and tem- 
perature data of version 3.3. Approximately 3500 H2O pro- 
files are retrieved on a daily basis from 316 hPa to 0.0002 hPa 
between 82° S and 82° N. MLS has a horizontal resolution of 
200-300 km along track, 7 km across track and 3-4 km in the 
vertical. In this study we will focus on the lower stratospheric 
pressure levels of lOOhPa, 82hPa, 68hPa, 56hPa, 46hPa, 
38hPa, 32 hPa, 26hPa, 22hPa and 18hPa on which MLS 
profiles are retrieved. The precision and accuracy of H2O at 
lOOhPa is 10% and 8%, respectively (Read et al., 2007). 
The MLS temperature product is retrieved from 316 hPa to 
0.001 hPa and interpolated on the same horizontal grid as 
H2O for the purpose of the study. It has a vertical resolution 
of about 5km at lOOhPa decreasing to 3km in the middle 
stratosphere. The precision at lOOhPa is 0.8 K whilst the un- 
certainty is about 2 K (warm bias to other instruments such 
as CHAMP and AIRS/AMSU Schwartz et al., 2008). MLS 

03 is measured using a radiometer at 240 GHz. Profiles are 
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reliable between the pressure levels 215 hPa and 0.02 hPa 
with a vertical resolution of 2.7 to 3 km from the upper tro- 
posphere to the mid-mesosphere (Froidevaux et al., 2008). 
To compare our observations with the QBO we use the Na- 
tional Centers for Environmental Prediction/National Center 
for Atmospheric Research (NCEP/NCAR) reanalysis daily 
zonal wind data, which are retrieved on 17 pressure levels 
from 1000 hPa to 10 hPa. The method used to determine ver- 
tical and meridional speeds of the BDC is similar to the one 
introduced by Schoeberl et al. (2008) and Flury et al. (2012). 
We calculate the time lagged correlation in the H 2 O time se- 
ries between two different levels. The average speed is then 
determined by the ratio of the distance between the two levels 
and the calculated optimal time lag. In order to obtain a time 
series of the average speed we use a one year time window 
for correlation and shift one level in time with respect to the 
other until a maximum correlation is reached. In choosing 
such a 1 yr time window for correlation seasonal variations 
cancel out and interannual variations are emphasised. This 
correlation procedure is repeated every day until the end of 
the measurement time series. Thus, each daily value of the 
derived speed time series represents the average speed over 
the following 365 days. For the vertical up welling speed we 
focus on the 68hPa and 56hPa levels at the Equator which 
are separated on the average by 1.12 km. The distance be- 
tween the two pressure levels is determined with the MLS 
version 3.3 geopotential height data (Schwartz et al., 2008), 
which were averaged over the respective time period of in- 
terest. We choose this altitude to make sure to be completely 
above tropical convection. To evaluate the speed of the shal- 
low meridional branch, we correlate the 100 hPa levels at the 
Equator and 40° latitude North and South, respectively. We 
use this speed as a proxy for the shallow branch knowing that 
the derived absolute speed is tracer dependent. We focus on 
the 100 hPa level because above 68 hPa H 2 O loses its distinc- 
tive tracer capabilities towards the mid-latitudes due to the so 
called tropical pipe (Plumb, 1996) which is stronger than in 
the lower stratosphere and inhibits mixing more effectively. 
We estimate the uncertainty of the derived speeds to be 5 %. 
The main contributor is the uncertainty in the determination 
of the optimal time lag, which is about 2 days in the here 
shown analysis. 

3 Results 

3.1 Water vapour transport 

Due to its long chemical lifetime of the order of tens of years 
in the lower stratosphere (Brasseur and Solomon, 2005) H 2 O 
can be used as a tracer for atmospheric transport given spatial 
gradients of concentrations in this region. Gradients exist in 
our region of interest as well in the vertical at the Equator 
as on the horizontal between tropics and mid-latitudes in the 
lower stratosphere. 


H 2 0 deviation from mean at 100.00 hPa 
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year 

Fig. 3. Relative deviation from the average zonal mean H 2 O as a 
function of time and latitude at lOOhPa. Red colours show above 
average concentrations, blue shows below average values. The vis- 
ible seasonal cycle with an amplitude of about 30 % is transported 
from the tropics towards the poles by the meridional BDC as indi- 
cated by the arrows and can be identified along the slant lines of 
same colours. The slope is steeper towards the NH which is indica- 
tive for a stronger and faster BDC. 


Tropical zonal mean H 2 O profile time series are suitable to 
visualise vertical transport mechanisms. Figure 2 shows the 
well known atmospheric tape recorder (Mote et al., 1996) 
as the pressure-time series of H 2 O deviations from its time 
mean at each pressure level. The seasonal cycle of H 2 O at 
lOOhPa is imprinted and transported upward by the ascend- 
ing branch of the BDC highlighted by the arrow. The signal 
is visible up to lOhPa. Short term interannual variations in 
H 2 O are produced by the QBO and El Nino Southern Oscil- 
lation (ENSO). Minimum values at 100 hPa occur during NH 
winter. The irregular interannual variability between 200 hPa 
and lOOhPa follows the ENSO signal. During El Nino the 
atmosphere is warmer and carries more H 2 O (red) e.g., in 
2005, 2007 and 2010. A strong La Nina occurred in 2007- 
2008 and 2010-2012 where H 2 O is well below the average 
(blue). The variable speed of the ascending air masses can be 
estimated by computing the slope of the alternating coloured 
surfaces above lOOhPa highlighted by the arrow. 

Figure 3 shows latitude-time series of H 2 O deviations 
from its zonal mean at lOOhPa. The seasonal cycle is clearly 
visible over the whole latitude range. Red colours represent 
higher than average H 2 O whilst blue shows lower than av- 
erage values. Furthermore, we observe that the alternating 
bands of positive (red) and negative (blue) H 2 O anomalies 
propagate towards the high-latitudes which is highlighted by 
the two arrows representing the shallow meridional branch of 
the BDC as previously shown in Fig. 1. The arrow is steeper 
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in the NH which means a faster transport, hence, a stronger 
BDC. The negative anomalies at high southern latitudes in 
the second half of each year stem from the very cold Antarc- 
tic polar vortex and the buildup of polar stratospheric clouds 
which dehydrate the 100 hPa level. 

3.2 BDC speed variations 

To compute the average speed of the vertical upwelling 
branch of the BDC as highlighted by the thick upward point- 
ing arrow in Fig. 1, we correlate the time series of the 68 hPa 
(~18.80km) and 56hPa (~19.92km) level in one year win- 
dows over the whole MLS time series. Correlation is strong 
with coefficients r ~ 0.9. Since the tape recorder signal is as- 
cending, the H 2 O at 56hPa is lagged in time with the one at 
the lower level. We determine the optimal time lag in order 
to get maximum correlation between the two levels. Thus, 
the ascent rate can be computed from the ratio of the mean 
distance between the two levels (A z ~ 1.12 km) and the op- 
timal time lag. The same calculation is done for the merid- 
ional component, which constitutes the shallow branch of the 
BDC. We compute the time lagged correlation between the 
100 hPa levels at the Equator and 40° N and 40° S (Ay ~ 
4440 km), respectively. Correlation coefficients are also gen- 
erally high with r ~ 0.9. The uncertainty of the method is 
composed of the uncertainty in the distance between the lev- 
els (~ 2 %) and the precision of the time lag (~ ±2 day). The 
two sum up to about 5 % uncertainty in the derived speed 
which is significantly smaller than the interannual variability 
of 21 % found in the study. As already pointed out in the in- 
troduction, we compute the speed of the poleward propagat- 
ing water vapour anomaly. Since this does not result from a 
pure transport process of mass, but includes two-way mixing 
the speed is tracer dependent. However, since poleward resid- 
ual mass transport is dominant (Trepte et al., 1993) our signal 
is a good proxy for the speed of the BDC and for simplicity 
we just call it speed of the meridional branch. Figure 4 shows 
the results of the mean upwelling- and meridional speeds of 
the BDC derived from the correlation of H 2 O time series. 
Note that the speed v in both hemispheres is taken positive 
towards either pole. The meridional transport in the NH (v, 

NH on the order of 1 ms -1 ) is about 5000 times faster than 
the vertical ascent (w) which is ~0.2mms _1 . The speed to- 
wards the SH (v, SH) is about half the one in the NH, which 
confirms the interpretation of Fig. 3 with a slower BDC in the 
SH. Also the austral variability of about 10 % is half the vari- 
ability in the NH (21 %). The variability in the vertical ascent 
(~ 21 %) is very similar to the one in the meridional to the 
NH. Furthermore, there is a prominent ~ 2 yr oscillation in 
unh and w which is linked to the well studied QBO. How- 
ever, the observed anti-correlation of the two BDC branches 
in the NH was rather surprising. When the vertical upwelling 
branch speeds up, the meridional slows down. Both are part 
of the BDC, but show opposite interannual variations. Sec- 
tion 3.4 will discuss these particular findings in more detail. 
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Fig. 4. Speed of the meridional branch of the BDC at lOOhPa to 
the NH (v, NH dashed black) and to the SH (v, SH, dashdotted 
black) indicated on the left axis. The equatorial vertical ascent rate 
(w, solid blue) at 68hPa is shown on the right axis. Speeds show 
a prominent ~ 2yr oscillation related to the QBO. Both v, NH and 
w show interannual variability of about 21 %. The vertical speed 
increases in phase with a meridional speed decrease and vice versa. 
The speed to the SH shows less variability and is out of phase with 
the NH for the most part, except for the first 3 yr where they are 
correlated. 


Since we use one year increments of H 2 O data for cor- 
relation we do not retrieve the typical seasonal variation of 
the Brewer-Dobson circulation which is picking up in the 
months September to March every year due to enhanced 
high-latitude stratospheric wave breaking. Niwano et al. 
(2003) showed time series of vertical ascent rate profiles us- 
ing HALOE H 2 O and CH 4 data where the seasonal cycle is 
visible with maxima occurring during NH winter. To verify 
the Niwano et al. (2003) results, we determine the time it 
takes for an air parcel to ascend by one pressure level using 
MLS daily zonal mean H 2 O data averaged between 8 ° S- 
8 ° N. We are specifically investigating the time it takes to 
shift a lower stratospheric segment of a profile up by one 
MLS pressure level as can be identified in Fig. 5. Thereby, 
we determine the number of days j it takes to obtain the best 
correlation of the lOOhPa to 56hPa profile at time ti (left 
blue) with a profile taken one MLS pressure level higher up 
i.e., from 82hPa to 46hPa at a time ti+j (right blue). This 
calculation is repeated for every day and subsequent pressure 
levels up to 10 hPa in order to get a time series of the average 
vertical ascent rate profile at a higher temporal resolution. 
We use MLS geopotential height data to compute the average 
height difference between the considered levels. The increas- 
ing spacing with altitude between the levels shown in Fig. 6 
(from Az = 1.11 km to Az = 1.31 km) is taken into account. 
The ascent rate is, thus, derived by the covered distance di- 
vided by the computed time lag. The speed is attributed to 
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Profile at time t Profile j days later at t 

i i+j 



Fig. 5. Schematics on how the vertical ascent rate is computed us- 
ing the Niwano et al. (2003) method. A portion of a H 2 O profile 
is taken at a time t[ between 100 hPa and 56hPa (left, blue). Due 
to BDC ascent the profile will be shifted upward some time later. 
At a specific time ti+j, j days later, the structure can be found one 
pressure level higher up (right, blue). In order to compute this time, 
we correlate the blue H 2 O segment on the left with the segment one 
level higher up (right, blue) on subsequent days until maximum cor- 
relation is reached after j days. The ratio of the ascended distance 
and the time lag j yields the average ascent rate. 


the centre pressure level of the considered altitude layer. Fig- 
ure 6 shows the time series of ascent rate profiles in the 
lower stratosphere between 68hPa (~ 18.8 km) and 18hPa 
(~ 27.2km). The speed varies between about 0.1 mm s -1 
and 0.5 mm s -1 with maximum values occurring usually dur- 
ing the NH autumn and winter months when the BDC is 
stronger. The values are greater than those obtained with 
the previous method because of the better time resolution 
which allows us to identify the seasonal cycle of the BDC. 
The average speed in the considered altitude range is about 
0.24 mm s -1 which again compares well to Fig. 4 taking into 
account the small increase of speed with altitude. However, 
to better compare the data with the vertical speed shown in 
Fig. 4, we compute the anomalies by removing the seasonal 
cycle which is shown in Fig. 7. Again, we observe a strong 
QBO signal which slowly descends with a speed of about 1- 
1.5 km month -1 from the upper boundary of the plot. Lower 
ascent rates are observed in 2006, 2008 and 2010-2011. The 
variability of about 50 % is bigger than observed in Fig. 4 and 
is due to the different method used to compute the speed. 

Table 1 shows an overview of our ascent rates in the lower 
stratosphere at six different pressure levels together with data 
published by Niwano et al. (2003), Mote et al. (1996) and 
Schoeberl et al. (2008). The values are extracted from the fig- 
ures in the corresponding publications. Our values represent 
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Fig. 6. Average vertical speed in the lower stratosphere between 
68hPa (~ 18.8 km) and 18hPa (~ 27.2 km) as a function of time 
calculated according to Niwano et al. (2003). A seasonal cycle is 
apparent with higher speeds (red) increasing with altitude and to- 
wards the end of each calendar year during NH winter, which is 
consistent with the characteristics of the BDC which is strongest 
from September to March. During NH summer the vertical speed 
decreases (blue) due to the reduced planetary wave breaking, which 
is the main driver of the BDC. 


Ascent Rate Anomaly (deseasoned) between 8°S-8°N % 
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Fig. 7. Anomalies of the ascent rates shown in Fig. 6 with the sea- 
sonal cycle removed. A distinct ~2 yr cycle is visible which is de- 
scending slowly at a rate of about 1.5 km month -1 . The variability 
is about 50 % with lower ascent rates (blue) observed in 2006, 2008 
and 2010-2011. The negative anomalies occur during the QBO 
westerly phase as will be shown in Fig. 10. Periods of faster up- 
welling are indicated in red. 
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Table 1 . Comparison of tropical ascent rates w [mm s -1 ] in the lower stratosphere from different studies. Niwano et al. (2003) and Mote et al. 
(1996) used HALOE satellite observations and averaged between ±12.5°, Schoeberl et al. (2008) used a combination of 15 yr of HALOE 
and Aura/MLS H 2 O observations to compute w at the Equator. They also displayed results of the GEOS-4 general circulation model. Ascent 
rates are increasing with altitude and very similar to each other. However, the GEOS-4 model underestimates the observations. The slow 
tropical ascent is in contrast to the faster meridional speed (~ 1ms -1 ). 


p [hPa] 

~ z [km] 

This Study 
w [mms -1 ] 

Niwano 03 
w [mms -1 ] 

Mote 96 
w [mms -1 ] 

Schoeberl 08 
w [mms -1 ] 

GEOS-4 GCM 
w [mms -1 ] 

70 

18.7 

0.20 



0.30 

0.12 

60 

19.6 

0.22 

0.22 

0.20 

0.30 

0.12 

50 

20.7 

0.22 

0.24 

0.22 

0.30 

0.15 

40 

22.0 

0.23 

0.25 

0.25 

0.32 

0.19 

30 

23.8 

0.27 

0.28 

0.30 

0.35 

0.25 

20 

26.4 

0.31 

0.33 

0.40 

0.50 

0.36 


Ascent rate and temperature anomaly at Equator and 68 hPa 



year 

Fig. 8. The vertical ascent rate at the Equator between 68 hPa and 
56 hPa (dashed black, left y-axis) is compared to the MLS tempera- 
ture anomaly (solid blue, right y-axis) at the same latitude and pres- 
sure. The temperature time series has been smoothed by a one year 
forward moving average to match the method which is used to de- 
termine the ascent rate. Each value represents the average over the 
next year. Temperature is anti-correlated with the vertical ascent as 
expected. 


the average of the data shown in Fig. 6 at corresponding pres- 
sure levels. All studies agree on the increase with altitude and 
the values agree well among each other. Taking into account 
the respective uncertainties (~ 5-25 %) given in the publi- 
cations the range of values overlap. Furthermore, Schoeberl 
et al. (2008) who obtained the highest values using a combi- 
nation of HALOE and Aura/MLS H 2 O show also the high- 
est standard deviation. In their work they compared the as- 
cent rates with the GEOS-4 general circulation model which 
shows overall the lowest values. The slow ascent rate is in 
contrast with the much faster meridional speed of ~ 1ms -1 . 


3.3 Link of temperature and vertical ascent 

The speed of the BDC circulation plays a key role in steering 
the tropical tropopause temperature as well as the tempera- 
ture in the lower stratosphere. An increase in speed of the 
BDC leads to a cooling of the tropical tropopause and the 
stratosphere as was already shown by Niwano et al. (2003); 
Randel et al. (2006); Yulaeva et al. (1994). Triggered by en- 
hanced up welling, the reaction of the temperature is a com- 
bination of adiabatic cooling and ozone decrease due to ad- 
vection of lower ozone concentrations from below (Yulaeva 
et al., 1994; Randel et al., 2006). Less ozone results in less 
radiative heating. In order to compare our time series of the 
ascent rate (see Fig. 4) with the temperature anomalies at the 
same latitude and pressure we need to compute a moving av- 
erage over a one year period since the calculated speeds rep- 
resent the average over the following year. Figure 8 shows the 
time series of the MLS temperature anomaly at the Equator 
and 68 hPa smoothed by a one year forward moving average 
in blue and the ascent rate in dashed black. Each point in the 
time series represents the average over the next year. As ex- 
pected, an increase of the vertical velocity is accompanied by 
a decrease in temperature. Clearly visible again is the QBO 
modulation producing the ~ 2 yr cycle. The expected anti- 
correlation between the temperature and the vertical ascent 
can be seen as a validation of our method to compute the 
speed of the BDC. 

3.4 QBO influence on BDC 

The QBO is associated with a secondary meridional circu- 
lation anomaly that does not only affect the residual trans- 
port to mid-latitudes, but also the vertical ascent in the trop- 
ical stratosphere (Reed, 1964; Plumb and Bell, 1982; Choi 
et al., 2002; Punge et al., 2009; Ribera et al., 2004). Fig- 
ure 9 shows schematically the QBO associated circulation 
anomaly as a latitude-height section on the left and the verti- 
cal profile of the equatorial zonal wind u on the right. Below 
the maximum easterlies (E) there is enhanced vertical ascent 


www.atmos-chem-phys.net/13/4563/2013/ 


Atmos. Chem. Phys., 13, 4563-4575, 2013 


4570 


T. Flury et al.: Variability of the Brewer-Dobson circulation 



Fig. 9. Sketch of the QBO secondary meridional circulation on the 
left and the zonal wind profile at the Equator on the right. Arrows in- 
dicate vertical and meridional motion with colours meaning warm- 
ing (red) and cooling (blue). Below the altitude of maximum easter- 
lies (E) air is rising faster in the easterly shear zone (ESZ). At the al- 
titude of maximum easterlies air is diverging towards the poles and 
increases the speed of the BDC, the opposite is found at the altitude 
of maximum westerlies (W) where air converges at the Equator (Eq) 
and leads to downdrafts below that level which decrease the speed 
of the ascent and warm the atmosphere inside the westerly shear 
zone (WSZ). From our results of the hemispheric asymmetry in the 
meridional BDC-branch we conclude that the secondary meridional 
circulation is less strong in the SH by a factor of 2 and reduced the 
arrows accordingly in size. Figure adapted from Choi et al. (2002); 
Plumb and Bell (1982); Punge et al. (2009). 


inside the easterly shear zone (EZS), which is associated 
with lower temperatures. At the altitude of maximum east- 
erlies (E) we find meridional divergence towards the poles. 
Inside the westerly shear zone (WSZ) just below the maxi- 
mum westerlies (W) we find a descent anomaly, which is as- 
sociated with higher temperatures and reduces the speed of 
the prevailing BDC ascent. At the height of maximum west- 
erlies we find a meridional convergence zone at the Equator. 
We adapted the sketch from Choi et al. (2002) and Punge 
et al. (2009) who drew two symmetrical cells separated by 
the Equator. However, Pena-Ortiz et al. (2008) report on an 
asymmetry in autumn and spring when there is only one cell 
which is shifted completely into the winter hemisphere. Our 
results also suggest an asymmetry in the meridional com- 
ponents since the interannual variability of the meridional 
branch towards the NH is twice as strong as to the SH (see 
Fig. 4). For this reason we adapted the size of the meridional 
arrows to highlight the asymmetry. To sum up, easterly winds 
are associated with cold anomalies and higher vertical ascent 
rates below the level of maximum easterlies whilst westerly 
winds are associated with warm anomalies and lower vertical 
ascent rates below the level of maximum westerlies. 

To compare the theory of Fig. 9 with MLS measurements 
we compare the ascent rate with the zonal wind of a level 
higher up. In our case we choose the ascent rate calculated 
between the levels 68hPa and 56hPa and the average zonal 


Equatorial zonal wind and vertical speed 



Fig. 10. Ascent rate w at the Equator between the 68 and 56 hPa 
level (dashed black, left y-axis) together with the average equatorial 
zonal wind u (solid blue, right y-axis) higher up at 32 hPa. In order 
to compare to the ascent rate the zonal wind is smoothed with a one 
year forward moving average. The data confirm the sketch in Fig. 9 
and show maximum ascent rates during the easterly phase ( u < 0) 
and minimum during the westerly phase (u >0). 

wind at the 32hPa level, both time series are taken at the 
Equator. The results are displayed in Fig. 10. The ascent rate 
is plotted in black (dashed) and varies between 0.14 mm s -1 
and 0.24 mm s _1 . The NCEP average zonal wind (solid blue) 
is anti-correlated (r = —0.84) with the ascent rate and con- 
firms the schematics of Fig. 9. The ascent rate starts increas- 
ing right after maximum westerlies- and decreases after the 
maximum easterlies are reached. The upward motion is slow- 
est during the westerly wind regime. Note that even during 
the westerlies the ascent rate is positive, the red downward 
arrow in Fig. 9 only indicates the QBO induced anomaly 
which slows down the prevailing upward motion but does 
not reverse it. 

Figure 1 1 shows the same zonal wind time series as the 
previous figure, but together with the average meridional 
speed between the Equator and 40° N at 100 hPa. The zonal 
wind speed correlates with the meridional speed, the correla- 
tion coefficient here is r = 0.73. This behaviour can also be 
explained with the sketch in Fig. 9. The considered 100 hPa 
level for meridional transport is well below the 32 hPa of the 
zonal wind and, for example, in the case of maximal eastward 
wind (E in Fig. 9) the converging arrows in the lowest part of 
the figure apply which slow down the transport towards the 
poles. During the westerly phase, the meridional transport 
gets faster below, which also applies to our case here. 
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Equatorial zonal wind and meridional speed to NH 



Fig. 11. Speed of the meridional branch (dashed black, left y-axis) 
of the BDC to the NH at lOOhPa and the average equatorial zonal 
wind at 32hPa (solid blue, right y-axis). The meridional speed is 
in phase with the wind and increases (decreases) as suggested by 
Fig. 9 during the westerly/positive (easterly/negative) phase of the 
zonal wind u. 

4 Discussion 

With our time lagged correlation method using water vapour 
time series we derived the speed of the vertical upwelling 
branch of the BDC at the Equator as well as a proxy for 
the speed of the shallow meridional branch. For reasons of 
simplicity, we use the term meridional speed for the derived 
signal although it is tied to the used tracer H 2 O. Whilst the 
tape recorder signal (tropical upwelling) represents a slow 
diabatic ascent and is due to vertical mass advection the 
meridional transport in the shallow branch is a combination 
of residual (poleward) mass transport and two-way mixing. 
The correlation of H 2 O at lOOhPa between 40° S/N and the 
Equator is, thus, more regarded as a proxy for the speed of 
the shallow meridional component rather than an absolute 
measure of the speed. The reason for this is two-way mixing 
and the considered tracer gradient field, which determines 
diffusion. However, since the considered propagation of the 
H 2 O anomaly between the Equator and 40° is predominantly 
due to poleward residual mass transport (Trepte et al., 1993), 
we consider it as a good proxy for the shallow branch of 
the BDC. Two-way mixing, which is faster in NH summer 
and autumn (Bonisch et al., 2009; Birner and Bonisch, 2011; 
Hegglin and Shepherd, 2007), can be observed in an O 3 time- 
latitude plot at lOOhPa shown in Fig. 12 where mid-latitude 
air with higher O 3 concentrations enters the tropics. Never- 
theless, during winter and spring the tropical lower strato- 
sphere is much more isolated from the extra-tropics and there 
is no visible mixing into the tropics. However, at the same 
time we clearly observe H 2 O transport towards the poles 
(Fig. 3). The mixing barrier can be identified at approx. ±20° 


Log(0 3 ) @ 100 hPa Log(0 3 [ppm]) 



year 


Fig. 12. Logarithm of zonal mean MLS O 3 volume mixing ratios 
(ppm) at lOOhPa. O 3 rich air from mid-latitudes (red-orange) is 
mixed into the tropics especially in NH summer and autumn. Dur- 
ing winter and spring the low O 3 concentrations (blue) in the tropics 
are more isolated because mixing equatorwards is strongly reduced. 


latitude and is referred to as the tropical pipe (Plumb, 1996). 
The pipe separates tropical- from mid-latitude air which can 
be seen in steep gradients of various trace gas distributions 
such as e.g., N 2 O, O 3 and H 2 O. The tropical pipe is stronger 
above 70 hPa than below where it is leaky especially during 
summer. 

In order to monitor transport processes using H 2 O as a 
tracer, we need an observable gradient, which exists on the 
lOOhPa level between the tropics and mid-latitudes. Mid- 
latitude H 2 O at 100 hPa originates in the tropics and is mostly 
transported through the shallow branch of the BDC. In late 
winter and spring measurements show an increase in the age 
of lower stratospheric air which is indicative for down welling 
from the deep branch of the BDC (Bonisch et al., 2009). 
However, this only happens over relatively short time peri- 
ods. The same authors used a concept of mass balance to 
compute the contribution of the deep and shallow branch to 
extratropical lower stratospheric air and demonstrated that 
downwelling plays a role only during spring whereas direct 
horizontal transport from the tropical tropopause layer con- 
tributes during the whole year. Since we use one year time 
segments of H 2 O to derive the meridional speed, the deep 
branch of the BDC does not influence our results signifi- 
cantly. Moreover, the meridional gradient between tropics 
and mid-latitudes is stronger than the vertical gradient around 
the lOOhPa level (100-82 hPa) so that the influence of mid- 
latitude vertical transport on our derived meridional speed is 
small. 

Our results show that the QBO has a significant influence 
on the vertical and meridional speed of the BDC. Below the 
zone of maximum easterlies (QBO east) the vertical speed 
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increases and the meridional speed decreases. The opposite 
happens below the zone of maximum westerlies (QBO west) 
as suggested by Fig. 9. Thus, the QBO has an influence on 
how the upwelling air masses are distributed into the trop- 
ical stratosphere. We suggest an influence on the tropical 
pipe (Plumb, 1996) to explain the variability of the speed to- 
wards mid-latitudes. Figure 13 sketches the tropical pipe and 
Brewer-Dobson transport for the QBO easterly phase on the 
left and the QBO westerly phase on the right. We assume 
a constant vertical supply of tropospheric air that crosses 
the tropical tropopause during both QBO phases. The en- 
hancement of the respective branches is highlighted by the 
increased and gray shaded arrows. During the easterly phase 
the meridional speed towards the NH is about 40 % smaller 
than during west. We suggest that during this time the trop- 
ical pipe is stronger and mixing with the mid-latitudes takes 
more time. On the other hand for continuity reasons the ver- 
tical ascent rate is increased by about 40 % (see Fig. 4) com- 
pared to the westerly phase. However, during the QBO west- 
erly phase we observe the opposite, meridional transport is 
faster and vertical speed decreases. We, thus, assume that the 
tropical pipe is more leaky and allows for faster transport 
towards the mid-latitudes which in turn reduces the vertical 
speed because the increased meridional flux leaves less mass 
to be moved upward. The reason could be that the lower al- 
titude limit of the tropical pipe varies as a function of the 
QBO. We assume that during QBO- west the tropical pipe 
starts higher up than during QB O-east which would explain 
the difference in the speed of mixing and transporting H 2 O 
towards mid-latitudes. 

The observed asymmetry between the meridional speed in 
the SH and NH is consistent with the theory of the QBO. Its 
influence is more visible in the NH because planetary wave 
activity at NH mid-latitudes, the primary driver of the BDC, 
depends strongly on the QBO because planetary wave prop- 
agation is a function of the variable zonal mean flow. Since 
planetary wave activity is lower in the SH the QBO influence 
on the meridional branch is small. However, as mentioned in 
the introduction, wave activity inside the tropical pipe (Ran- 
del et al., 2008; Kerr-Munslow and Norton, 2006; Norton, 
2006; Garny et al., 2011; Engida and Folkins, 2012) has also 
an influence on tropical upwelling and the relative impor- 
tance of both mechanisms is still a matter of debate. We sug- 
gest that modelling groups try to reproduce our observations 
( w and v) with the goal to quantify the relative contribution 
of tropical and high-latitude wave activity on the BDC by 
either switching on or off the different wave forcings. 

Our results show that the QBO is the main driver of in- 
terannual variability in upwelling and meridional transport 
whilst the BDC controls the seasonal variability. On the other 
hand, the influence of the troposphere on the stratosphere 
is dominated by ENSO, which modulates the temperature 
and water vapour content of the troposphere particularly the 
tropical tropopause temperature. Nevertheless, we do not dis- 
cuss ENSO in more detail since its influence on interannual 
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Fig. 13. Interpretation of the QBO influence on the BDC variability. 
During the QBO easterly phase (left) the tropical pipe is stronger 
and it takes more time for air masses to cross the barrier towards 
mid-latitudes. Thus, the meridional speed is slower. For continu- 
ity, the vertical branch has to speed up which is shown by the big 
shaded vertical arrow. On the other hand, during the QBO westerly 
phase the tropical pipe is shifted up, more leaky and transport to- 
wards the mid-latitudes faster (right). Again, for continuity reasons, 
the vertical transport has to slow down which agrees well with our 
results. 

variability of the BDC is weaker than the QBO’s. However, 
we note that ENSO and QBO get out of phase in 2009/2010 
(Liang et al., 2011) which corresponds also to the time when 
the vertical ( w ) and meridional (vNH) velocity lose their an- 
ticorrelation. Whenever ENSO and QBO occur in phase the 
signatures on TTL- and lower stratospheric water vapour are 
the most prominent. Yet, a detailed analysis would require to 
understand the effects of the interaction of QBO and ENSO 
which is beyond the scope of this study. 


5 Conclusions 

We used the water vapour tape recorder signal to derive the 
speed of the tropical upwelling branch of the BDC from 7 yr 
of Aura/MLS measurements. In a similar way, we computed 
a proxy for the speed of the shallow meridional branch of 
the BDC in the lowermost stratosphere of both hemispheres. 
We found that the interannual variability of the BDC speed 
depends on the phase of the QBO. The tropical upwelling 
branch as well as the meridional branch towards the NH show 
interannual amplitudes of about 21 % whereas the SH merid- 
ional branch is less variable. The average speed to the NH 
is 1.15 ms -1 with an amplitude of 0.24ms -1 whilst it is 
only 0.58 ms -1 with an amplitude of 0.06 ms -1 to the SH. 
The tape recorder signal ascends with an average speed of 
0.20 mm s -1 with an amplitude of 0.04 mm s -1 in the lower 
stratosphere at the Equator. When using the method to re- 
trieve short term variability of the vertical ascent we found 
maximum values of 0.5 mm s -1 during NH winter higher 
up at 20 hPa. Moreover, we found that the variability in the 
meridional (NH) and vertical branch are anti-correlated. This 
anti-correlation is regulated by the so called QBO secondary 
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meridional circulation which is sketched in Fig. 9. We sug- 
gest that during the QBO easterly phase the tropical pipe is 
stronger and it takes more time for tropical air masses to 
mix into mid-latitude air. However, during the QBO west- 
erly phase the tropical pipe appears more leaky and merid- 
ional transport is faster. We suggest that the lower altitude 
limit of the tropical pipe is higher up during QBO- west. Our 
method was able to show the asymmetry of NH and SH and 
the different influence of the QBO on the BDC. If used by 
modellers our results may help to quantify the relative im- 
portance of tropical wave dissipation and high-latitude wave 
breaking on the BDC. 

MLS temperature measurements confirm the anti- 
correlation between temperatures and vertical speed and 
show cooling (warming) during the easterly (westerly) 
shear phase while vertical speed is higher (lower) and 
meridional speed is lower (higher) (Fig. 8). NCEP/NCAR 
zonal wind reanalysis data also show the influence of the 
QBO on the tropical up welling and meridional branch of 
the BDC. Strong easterlies lead to faster vertical ascent 
whilst westerlies reduce the ascent rate (Fig. 10). On the 
other hand, westerlies higher up increase the meridional 
speed at the 100 hPa level, whereas it is lower during the 
easterly phase (Fig. 11). However, since the QBO secondary 
meridional circulation is a function of altitude and time 
the conclusions derived here are based largely on MLS 
observations at lOOhPa and 68hPa. We suggest that further 
studies concentrate on the meridional branch of the BDC 
higher up in the stratosphere using H2O and methane (CH4) 
data. H2O alone is less suitable as a tracer for horizontal 
transport above 68 hPa due to the oxidation of CH4 which 
is a significant source of H2O in the stratosphere. General 
circulation models can be tested for the correct representa- 
tion of the atmospheric tape recorder. The underlying H2O 
distribution in the upper troposphere and the seasonal cycle 
is crucial for a correct representation, a task where most 
models still perform poorly as noted by Jiang et al. (2012). 
It is also difficult to find the QBO signal in the modelled 
tape recorder (J. Jiang, personal communication, 2012). 
MLS has been a very stable instrument for monitoring 
middle atmospheric dynamics though its primary strength 
is on atmospheric composition and chemistry. It is, thus, 
highly recommended to be used intensively for validation 
and improvements of global chemistry climate models since 
important atmospheric oscillations such as ENSO and QBO 
are well recorded by MLS measurements. 
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